Evidence-based Decision Making: Infectious Disease Modeling Training for Policymakers in East Africa.

Background: Mathematical modeling of infectious diseases is an important decision-making tool for outbreak control. However, in Africa, limited expertise reduces the use and impact of these tools on policy. Therefore, there is a need to build capacity in Africa for the use of mathematical modeling to inform policy. Here we describe our experience implementing a mathematical modeling training program for public health professionals in East Africa. Methods: We used a deliverable-driven and learning-by-doing model to introduce trainees to the mathematical modeling of infectious diseases. The training comprised two two-week in-person sessions and a practicum where trainees received intensive mentorship. Trainees evaluated the content and structure of the course at the end of each week, and this feedback informed the strategy for subsequent weeks. Findings: Out of 875 applications from 38 countries, we selected ten trainees from three countries - Rwanda (6), Kenya (2), and Uganda (2) - with guidance from an advisory committee. Nine trainees were based at government institutions and one at an academic organization. Participants gained skills in developing models to answer questions of interest and critically appraising modeling studies. At the end of the training, trainees prepared policy briefs summarizing their modeling study findings. These were presented at a dissemination event to policymakers, researchers, and program managers. All trainees indicated they would recommend the course to colleagues and rated the quality of the training with a median score of 9/10. Conclusions: Mathematical modeling training programs for public health professionals in Africa can be an effective tool for research capacity building and policy support to mitigate infectious disease burden and forecast resources. Overall, the course was successful, owing to a combination of factors, including institutional support, trainees' commitment, intensive mentorship, a diverse trainee pool, and regular evaluations.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

A standardised differential privacy framework for epidemiological modeling with mobile phone data.

During the COVID-19 pandemic, the use of mobile phone data for monitoring human mobility patterns has become increasingly common, both to study the impact of travel restrictions on population movement and epidemiological modeling. Despite the importance of these data, the use of location information to guide public policy can raise issues of privacy and ethical use. Studies have shown that simple aggregation does not protect the privacy of an individual, and there are no universal standards for aggregation that guarantee anonymity. Newer methods, such as differential privacy, can provide statistically verifiable protection against identifiability but have been largely untested as inputs for compartment models used in infectious disease epidemiology. Our study examines the application of differential privacy as an anonymisation tool in epidemiological models, studying the impact of adding quantifiable statistical noise to mobile phone-based location data on the bias of ten common epidemiological metrics. We find that many epidemiological metrics are preserved and remain close to their non-private values when the true noise state is less than 20, in a count transition matrix, which corresponds to a privacy-less parameter ϵ = 0.05 per release. We show that differential privacy offers a robust approach to preserving individual privacy in mobility data while providing useful population-level insights for public health. Importantly, we have built a modular software pipeline to facilitate the replication and expansion of our framework.

Measurably recombining malaria parasites.

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis.

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Human Trypanosoma cruzi chronic infection leads to individual level steady-state parasitemia: Implications for drug-trial optimization in Chagas disease.

Currently available drugs against Trypanosoma cruzi infection, which causes 12000 deaths annually, have limitations in their efficacy, safety and tolerability. The evaluation of therapeutic responses to available and new compounds is based on parasite detection in the bloodstream but remains challenging because a substantial proportion of infected individuals have undetectable parasitemia even when using diagnostic tools with the highest accuracy. We characterize parasite dynamics which might impact drug efficacy assessments in chronic Chagas by analyzing pre- and post-treatment quantitative-PCR data obtained from blood samples collected regularly over a year. We show that parasitemia remains at a steady-state independently of the diagnostic sensitivity. This steady-state can be probabilistically quantified and robustly predicted at an individual level. Furthermore, individuals can be assigned to categories with distinct parasitological status, allowing a more detailed evaluation of the efficacy outcomes and adjustment for potential biases. Our analysis improves understanding of parasite dynamics and provides a novel background for optimizing future drug efficacy trials in Chagas disease. Trial Registration: original trial registered with ClinicalTrials.gov, number NCT01489228.

Declines in prevalence alter the optimal level of sexual investment for the malaria parasite Plasmodium falciparum.

Successful infectious disease interventions can result in large reductions in parasite prevalence. Such demographic change has fitness implications for individual parasites and may shift the parasite's optimal life history strategy. Here, we explore whether declining infection rates can alter Plasmodium falciparum's investment in sexual versus asexual growth. Using a multiscale mathematical model, we demonstrate how the proportion of polyclonal infections, which decreases as parasite prevalence declines, affects the optimal sexual development strategy: Within-host competition in multiclone infections favors a greater investment in asexual growth whereas single-clone infections benefit from higher conversion to sexual forms. At the same time, drug treatment also imposes selection pressure on sexual development by shortening infection length and reducing within-host competition. We assess these models using 148 P. falciparum parasite genomes sampled in French Guiana over an 18-y period of intensive intervention (1998 to 2015). During this time frame, multiple public health measures, including the introduction of new drugs and expanded rapid diagnostic testing, were implemented, reducing P. falciparum malaria cases by an order of magnitude. Consistent with this prevalence decline, we see an increase in the relatedness among parasites, but no single clonal background grew to dominate the population. Analyzing individual allele frequency trajectories, we identify genes that likely experienced selective sweeps. Supporting our model predictions, genes showing the strongest signatures of selection include transcription factors involved in the development of P. falciparum's sexual gametocyte form. These results highlight how public health interventions impose wide-ranging selection pressures that affect basic parasite life history traits.

Connecting the dots: understanding how human mobility shapes TB epidemics.

Tuberculosis (TB) remains a leading infectious cause of death worldwide. Reducing TB infections and TB-related deaths rests ultimately on stopping forward transmission from infectious to susceptible individuals. Critical to this effort is understanding how human host mobility shapes the transmission and dispersal of new or existing strains of Mycobacterium tuberculosis (Mtb). Important questions remain unanswered. What kinds of mobility, over what temporal and spatial scales, facilitate TB transmission? How do human mobility patterns influence the dispersal of novel Mtb strains, including emergent drug-resistant strains? This review summarizes the current state of knowledge on mobility and TB epidemic dynamics, using examples from three topic areas, including inference of genetic and spatial clustering of infections, delineating source-sink dynamics, and mapping the dispersal of novel TB strains, to examine scientific questions and methodological issues within this topic. We also review new data sources for measuring human mobility, including mobile phone-associated movement data, and discuss important limitations on their use in TB epidemiology.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: A malaria case study.

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for the genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, such as by measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ("4CAST") composed of highly diverse antigens, and a 129-locus panel ("AMPLseq") composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitaemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for the genomic epidemiology of infectious disease.

Evaluating the reliability of mobility metrics from aggregated mobile phone data as proxies for SARS-CoV-2 transmission in the USA: a population-based study.

BACKGROUND: In early 2020, the response to the SARS-CoV-2 pandemic focused on non-pharmaceutical interventions, some of which aimed to reduce transmission by changing mixing patterns between people. Aggregated location data from mobile phones are an important source of real-time information about human mobility on a population level, but the degree to which these mobility metrics capture the relevant contact patterns of individuals at risk of transmitting SARS-CoV-2 is not clear. In this study we describe changes in the relationship between mobile phone data and SARS-CoV-2 transmission in the USA. METHODS: In this population-based study, we collected epidemiological data on COVID-19 cases and deaths, as well as human mobility metrics collated by advertisement technology that was derived from global positioning systems, from 1396 counties across the USA that had at least 100 laboratory-confirmed cases of COVID-19. We grouped these counties into six ordinal categories, defined by the National Center for Health Statistics (NCHS) and graded from urban to rural, and quantified the changes in COVID-19 transmission using estimates of the effective reproduction number (Rt) between Jan 22 and July 9, 2020, to investigate the relationship between aggregated mobility metrics and epidemic trajectory. For each county, we model the time series of Rt values with mobility proxies. FINDINGS: We show that the reproduction number is most strongly associated with mobility proxies for change in the travel into counties (0·757 [95% CI 0·689 to 0·857]), but this relationship primarily holds for counties in the three most urban categories as defined by the NCHS. This relationship weakens considerably after the initial 15 weeks of the epidemic (0·442 [-0·492 to -0·392]), consistent with the emergence of more complex local policies and behaviours, including masking. INTERPRETATION: Our study shows that the integration of mobility metrics into retrospective modelling efforts can be useful in identifying links between these metrics and Rt. Importantly, we highlight potential issues in the data generation process for transmission indicators derived from mobile phone data, representativeness, and equity of access, which must be addressed to improve the interpretability of these data in public health. FUNDING: There was no funding source for this study.

Distinguishing gene flow between malaria parasite populations.

Measuring gene flow between malaria parasite populations in different geographic locations can provide strategic information for malaria control interventions. Multiple important questions pertaining to the design of such studies remain unanswered, limiting efforts to operationalize genomic surveillance tools for routine public health use. This report examines the use of population-level summaries of genetic divergence (FST) and relatedness (identity-by-descent) to distinguish levels of gene flow between malaria populations, focused on field-relevant questions about data size, sampling, and interpretability of observations from genomic surveillance studies. To do this, we use P. falciparum whole genome sequence data and simulated sequence data approximating malaria populations evolving under different current and historical epidemiological conditions. We employ mobile-phone associated mobility data to estimate parasite migration rates over different spatial scales and use this to inform our analysis. This analysis underscores the complementary nature of divergence- and relatedness-based metrics for distinguishing gene flow over different temporal and spatial scales and characterizes the data requirements for using these metrics in different contexts. Our results have implications for the design and implementation of malaria genomic surveillance studies.

Low parasite connectivity among three malaria hotspots in Thailand.

Identifying sources and sinks of malaria transmission is critical for designing effective intervention strategies particularly as countries approach elimination. The number of malaria cases in Thailand decreased 90% between 2012 and 2020, yet elimination has remained a major public health challenge with persistent transmission foci and ongoing importation. There are three main hotspots of malaria transmission in Thailand: Ubon Ratchathani and Sisaket in the Northeast; Tak in the West; and Yala in the South. However, the degree to which these hotspots are connected via travel and importation has not been well characterized. Here, we develop a metapopulation model parameterized by mobile phone call detail record data to estimate parasite flow among these regions. We show that parasite connectivity among these regions was limited, and that each of these provinces independently drove the malaria transmission in nearby provinces. Overall, our results suggest that due to the low probability of domestic importation between the transmission hotspots, control and elimination strategies can be considered separately for each region.

A nowcasting framework for correcting for reporting delays in malaria surveillance.

Time lags in reporting to national surveillance systems represent a major barrier for the control of infectious diseases, preventing timely decision making and resource allocation. This issue is particularly acute for infectious diseases like malaria, which often impact rural and remote communities the hardest. In Guyana, a country located in South America, poor connectivity among remote malaria-endemic regions hampers surveillance efforts, making reporting delays a key challenge for elimination. Here, we analyze 13 years of malaria surveillance data, identifying key correlates of time lags between clinical cases occurring and being added to the central data system. We develop nowcasting methods that use historical patterns of reporting delays to estimate occurred-but-not-reported monthly malaria cases. To assess their performance, we implemented them retrospectively, using only information that would have been available at the time of estimation, and found that they substantially enhanced the estimates of malaria cases. Specifically, we found that the best performing models achieved up to two-fold improvements in accuracy (or error reduction) over known cases in selected regions. Our approach provides a simple, generalizable tool to improve malaria surveillance in endemic countries and is currently being implemented to help guide existing resource allocation and elimination efforts.

The association between gold mining and malaria in Guyana: a statistical inference and time-series analysis.

BACKGROUND: Guyana reported a significant rise in malaria between 2008 and 2014. As there was no evidence of impairment of national malaria control strategies, public health authorities attributed the surge to a temporal increase in gold mining activity in forested regions. However, systematic analysis of this association is lacking because of the difficulties associated with collecting reliable data for both malaria and mining. We aimed to investigate the association between the international gold price and Plasmodium falciparum malaria transmission in Guyana between 2007 and 2019. We also aimed to evaluate the association between P falciparum cases and the El Niño-Southern Oscillation pattern, which has previously been suggested as a major driver of malaria. METHODS: We used national malaria surveillance data from Guyana to estimate the correlation over time between the international gold price and reported P falciparum infections in individuals who were likely to be involved in mining activities (ie, men and boys aged between 15 and 50 years who were living in mining regions) for each month between 2007 and 2019. We compared the estimates with those obtained from individuals who were unlikely to be directly involved in mining activities (ie, women, children aged 12 years and younger, and adults aged over 70 years) and estimates obtained from individuals living in non-mining regions. We also evaluated the correlation between P falciparum infections and the El Niño-Southern Oscillation pattern in the same subpopulations and time period. Lastly, we evaluated the performance of a statistical model formulated to estimate P falciparum infections in real time using the international gold price as the predictor variable. FINDINGS: The proportion of P falciparum malaria cases in temporary residents, which was used as a proxy for circulating individuals involved in gold mining, was highest during the years of peak gold price (ie, between 2008 and 2014). Cases of malaria in all demographic groups showed a strong positive correlation with the gold price, but only in regions with mining camps (0·88 [95% CI 0·84-0·89] for boys and men aged between 15 and 50 years and 0·80 [0·73-0·85] for the aggregated population of women, children aged 12 years and younger, and adults older than 70 years). The highest correlation occurred earlier in men and boys aged between 15 and 50 years, the demographic most likely to be miners, suggesting that transmission in mining camps is followed by infections in the community. On the basis of these findings, we were able to reliably forecast P falciparum malaria trends using only the gold price as the predictor variable. A 1% increase in gold price was associated with a 2·13% increase in P falciparum infections after 1 month in the mining populations, and with a 1·63% increase after 2 months in the non-mining populations. Lastly, La Niña climatic events showed an additional, smaller positive correlation with malaria transmission. INTERPRETATION: Our analysis provides evidence that the P falciparum malaria surge observed in Guyana between 2008 and 2014 was likely to have been driven mainly by an increase in gold mining, while climate factors might have contributed synergistically. We propose that the international gold price over time is a useful indicator of malaria trends. We conclude that the feasibility of malaria elimination in Guyana, and in other areas in the Amazon where malaria and gold mining overlap, should be evaluated against the challenges posed by rapidly rising gold prices. FUNDING: Ramón Areces Foundation, National Institutes of Health, and National Institute of General Medical Sciences.

Characterizing human mobility patterns in rural settings of sub-Saharan Africa.

Human mobility is a core component of human behavior and its quantification is critical for understanding its impact on infectious disease transmission, traffic forecasting, access to resources and care, intervention strategies, and migratory flows. When mobility data are limited, spatial interaction models have been widely used to estimate human travel, but have not been extensively validated in low- and middle-income settings. Geographic, sociodemographic, and infrastructure differences may impact the ability for models to capture these patterns, particularly in rural settings. Here, we analyzed mobility patterns inferred from mobile phone data in four Sub-Saharan African countries to investigate the ability for variants on gravity and radiation models to estimate travel. Adjusting the gravity model such that parameters were fit to different trip types, including travel between more or less populated areas and/or different regions, improved model fit in all four countries. This suggests that alternative models may be more useful in these settings and better able to capture the range of mobility patterns observed.

Genomics, social media and mobile phone data enable mapping of SARS-CoV-2 lineages to inform health policy in Bangladesh.

Genomics, combined with population mobility data, used to map importation and spatial spread of SARS-CoV-2 in high-income countries has enabled the implementation of local control measures. Here, to track the spread of SARS-CoV-2 lineages in Bangladesh at the national level, we analysed outbreak trajectory and variant emergence using genomics, Facebook 'Data for Good' and data from three mobile phone operators. We sequenced the complete genomes of 67 SARS-CoV-2 samples (collected by the IEDCR in Bangladesh between March and July 2020) and combined these data with 324 publicly available Global Initiative on Sharing All Influenza Data (GISAID) SARS-CoV-2 genomes from Bangladesh at that time. We found that most (85%) of the sequenced isolates were Pango lineage B.1.1.25 (58%), B.1.1 (19%) or B.1.36 (8%) in early-mid 2020. Bayesian time-scaled phylogenetic analysis predicted that SARS-CoV-2 first emerged during mid-February in Bangladesh, from abroad, with the first case of coronavirus disease 2019 (COVID-19) reported on 8 March 2020. At the end of March 2020, three discrete lineages expanded and spread clonally across Bangladesh. The shifting pattern of viral diversity in Bangladesh, combined with the mobility data, revealed that the mass migration of people from cities to rural areas at the end of March, followed by frequent travel between Dhaka (the capital of Bangladesh) and the rest of the country, disseminated three dominant viral lineages. Further analysis of an additional 85 genomes (November 2020 to April 2021) found that importation of variant of concern Beta (B.1.351) had occurred and that Beta had become dominant in Dhaka. Our interpretation that population mobility out of Dhaka, and travel from urban hotspots to rural areas, disseminated lineages in Bangladesh in the first wave continues to inform government policies to control national case numbers by limiting within-country travel.

Sociodemographic characteristics of missing data in digital phenotyping.

The ubiquity of smartphones, with their increasingly sophisticated array of sensors, presents an unprecedented opportunity for researchers to collect longitudinal, diverse, temporally-dense data about human behavior while minimizing participant burden. Researchers increasingly make use of smartphones for "digital phenotyping," the collection and analysis of raw phone sensor and log data to study the lived experiences of subjects in their natural environments using their own devices. While digital phenotyping has shown promise in fields such as psychiatry and neuroscience, there are fundamental gaps in our knowledge about data collection and non-collection (i.e., missing data) in smartphone-based digital phenotyping. In this meta-study using individual-level data from six different studies, we examined accelerometer and GPS sensor data of 211 participants, amounting to 29,500 person-days of observation, using Bayesian hierarchical negative binomial regression with study- and user-level random intercepts. Sensitivity analyses including alternative model specification and stratified models were conducted. We found that iOS users had lower GPS non-collection than Android users. For GPS data, rates of non-collection did not differ by race/ethnicity, education, age, or gender. For accelerometer data, Black participants had higher rates of non-collection, but rates did not differ by sex, education, or age. For both sensors, non-collection increased by 0.5% to 0.9% per week. These results demonstrate the feasibility of using smartphone-based digital phenotyping across diverse populations, for extended periods of time, and within diverse cohorts. As smartphones become increasingly embedded in everyday life, the insights of this study will help guide the design, planning, and analysis of digital phenotyping studies.

Deconstructing the parasite multiplication rate of Plasmodium falciparum.

Epidemiological indicators describing population-level malaria transmission dynamics are widely used to guide policy recommendations. However, the determinants of malaria outcomes within individuals are still poorly understood. This conceptual gap partly reflects the fact that there are few indicators that robustly predict the trajectory of individual infections or clinical outcomes. The parasite multiplication rate (PMR) is a widely used indicator for the Plasmodium intraerythrocytic development cycle (IDC), for example, but its relationship to clinical outcomes is complex. Here, we review its calculation and use in P. falciparum malaria research, as well as the parasite and host factors that impact it. We also provide examples of metrics that can help to link within-host dynamics to malaria clinical outcomes when used alongside the PMR.